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ABSTRACT 

We investigate the potential of submm-mm and submm-mm-radio photometric red- 
shifts using a sample of mm-selected sources as seen at 250, 350 and 500 /im by the 
SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources 
with reliable radio identifications in the GOODS-N and Lockman Hole North fields 
46 (73 per cent) are found to have detections in at least one SPIRE band. We explore 
the observed submm/mm colour evolution with redshift, finding that the colours of 
mm-sources are adequately described by a modified blackbody with constant optical 
depth r = {v/vqY where /3 = +1.8 and = c/100 /mi. We find a tight correla- 
tion between dust temperature and IR luminosity. Using a single model of the dust 
temperature and IR luminosity relation we derive photometric redshift estimates for 
the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spec- 
troscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm 
photometric redshifts offer a redshift accuracy of | Ax|/ (1+z) = 0.16 (< | Az\ >= 0.51). 
Including constraints from the radio-far IR correlation the accuracy is improved to 
|Az|/(l + z) = 0.15 (< |Az| >= 0.45). We estimate the redshift distribution of mm- 
selected sources finding a significant excess at z > 3 when compared to ~ 850 /tm 
selected samples. 

Key words: 


1 INTRODUCTION 

As is the case in the local Universe, star formation at high 
redshift takes place in regions rich in molecular gas and dust. 
Franceschini et al. (1991) and Blain & Longair (1993) thus 
anticipated that the distant galaxies hosting the majority of 
this activity would be discovered in the submm/mm wave- 
band, which benefits from a strong, negative A-correction 
for the rest-frame far-infrared (-IR) portion of their spectral 
energy distributions (SEDs). Observing at 450 and 850 /tm 
the SCUBA submm camera on the JCMT (Holland et al. 
1999) made it possible to detect the first examples of this 
galaxy population and represented a fundamental turning 
point in our understanding of galaxy formation in the dis- 
tant Universe (Smail, Ivison & Blain 1997; Barger et al. 
1998; Hughes et al. 1998; Eales et al. 1999; Dey et al. 1999). 

Early work on the submm galaxy population was of- 
ten hindered by small numbers of sources, and significant 
difficulties in finding multiwavelength counterparts given 
the relatively low angular resolution (FWHM> 10 arcsec) 
of ground-based submm/mm facilities. This situation has 
been much improved by a large increase in the number 
of submm/mm facilities, larger scale submm surveys (e.g. 
SHADES, Mortier et al. 2005; LESS, Weifi et al. 2009), im- 
provements in the quality of overlapping multiwavelength 
data (in particular radio and mid-IR imaging), as well as 
advances in cross-matching algorithms (Ivison et al. 2007; 
Chapin et al. 2009; Roseboom et al. 2009; Chapin et al. 
2011). As a result, our understanding of these galaxies has 
accelerated in recent years (e.g. Chapman et al. 2005; Greve 
et al. 2005; Alexander et al. 2005; Swinbank et al. 2008; Tac- 
coni et al. 2008; Pope et al. 2008; Menendez-Delmestre et 
al. 2009). 

The advent of the SPIRE instrument (Griffin et al. 
2010) on board Hersche 0 (Pilbratt et al. 2010) promises to 
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1 Herschel is an ESA space observatory with science instruments 
provided by Principal Investigator consortia. It is open for pro- 


shed new light on many issues relating to submm galaxies. In 
photometry mode, SPIRE operates in three submm bands 
centred on 250, 350 and 500 /tm, and thus the large scale 
surveys conducted with SPIRE will detect many 10,000s 
of submm galaxies (e.g. Eales et al. 2010; Oliver et al., in 
prep). SPIRE observations of existing submm/mm survey 
fields provides the crucial SED information shortward of 
500 /tm which cannot be routinely accessed from the ground 
(c.f. Kovacs et al. 2006), allowing us to characterise their IR 
luminosity (Lir) and dust temperature (7b) with unprece- 
dented accuracy (e.g. Chapman et al. 2010; Magnelli ct al. 
2010). 

Armed with this new information we can re-evaluate the 
potential of using submm (and mm) wavelength data alone 
to estimate the redshifts of submm galaxies. As submm 
galaxies tend to be both high redshift (< 2 2; Chapman 

et al. 2005) and highly dust obscured they are extremely 
faint at optical/near-IR wavelengths and hence obtaining 
redshifts via optical spectroscopy is challenging. Because 
of this the prospects for submm/mm photometric redshifts 
have been discussed many times since the early SCUBA 
surveys (Blain 1999; Hughes et al. 2002; Aretxaga et al. 
2003; Blain, Barnard & Chapman 2003; Pope & Chary 2010; 
Schulz et al. 2010), with mixed results. The main obstacle for 
submm/mm photometric redshifts is that redshift and dust 
temperature are degenerate for the simplest, single temper- 
ature, modified blackbody SED. To break this degeneracy 
some other constraint is needed, such as the IR luminosity 
to dust temperature relation (Blain 1999; Blain, Barnard & 
Chapman 2003) or information at other wavelengths, phys- 
ically related to the far IR emission (e.g. radio or mid-IR; 
Carilli & Yun 1999; Aretxaga et al. 2005; Aretxaga et al. 
2007). 

In this paper we investigate the potential of submm/mm 
photometric redshifts using a sample selected at SPIRE and 


posals for observing time from the worldwide astronomical com- 
munity. 
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mm-wavelengths. We focus on two key survey fields; Great 
Observatories Origins Deep Survey North (GOODS-N) and 
Lockman Hole North, which contain some of the deepest 
imaging on the sky at mm, submm, and radio wavelengths. 
In particular the Owen et al. (2008) and Morrison et al. 
(2010) VLA imaging of the Lockman Hole, and GOODS-N 
fields, represent the deepest ‘blank-held’ images currently 
available at 1.4 GHz. The depth of the radio imaging is cru- 
cial as without the positional information from identifica- 
tions in high angular resolution radio imaging it is often 
not possible to correctly deblend the highly confused SPIRE 
imaging. 

In ; [2] we introduce the parent mm-wavelength cata- 
logues utilised in this work. Sj3] describes the process under- 
taken to find multiwavelength identifications for the mm and 
SPIRE sources, U5] presents the SPIRE detection statistics 
for our mm-wavelength sample, and describes the proper- 
ties of those sources with both mm and SPIRE detections. 
Finally, [J5] considers the scientific implications of these re- 
sults, and we present our conclusions in Sj6] Throughout we 
assume a ACDM cosmology with ft a = 0.7, = 0.3 and 

No = 70 kins -1 Mpc -1 . 


2 DATA 

The starting point for this study is the available mm- 
wavelength imaging in the GOODS-N and Lockman Hole 
North field. 

In GOODS-N we make use of the combined MAMBO- 
AzTEC image and source catalogue of Penner et al. (2011). 
Combining the pre-existing 1.1-mm AzTEC (Perera et al. 
2008) and MAMBO 1.2-mm (Greve et al. 2008) imag- 
ing of GOODS-N the resulting image has a typical noise 
of 0.5 mjy at an effective wavelength of 1.16 mm over an 
area of 0.08 deg 2 . This combined image yields a list of 41 
1.16-mm sources detected at a significance of > 3.8<r and 
Si. i6 mm =2 — 10 mjy. In Lockman Hole North we make use 
of the recent MAMBO 1.2-mm imaging and source catalogue 
of Lindner et al. (2011). Lindner et al. present a list of 41 1.2- 
mm sources detected with S/N> 4 and Si. 2 mm = 2 — 5 mjy. 

In addition we make use of the SPIRE imaging at 250, 
350 and 500 pm obtained as part of the Science Demonstra- 
tion Phase (SDP) of Herschel by the Herschel Multi-Tiered 
Extragalactic Survey (HerMEqJ Oliver et al., in prep.). 
The SPIRE instrument, its in-orbit performance and its sci- 
entific capabilities are described by Griffin et al. (2010); its 
calibration methods and accuracy are outlined in Swinyard 
et al. (2010). We use images produced as described in Lev- 
enson et al. (2010). The SPIRE instrument has a FWHM of 
18, 25 and 36 arcsec at 250, 350 and 500 pm, respectively. 
HerMES imaging of GOODS-N and Lockman Hole North 
has a typical 5-cr depth at all SPIRE wavelengths of ~5 mjy 
and ~10mJy, respectively, ignoring the dominant contribu- 
tion from confusion (Smith et al., in press). 

At radio wavelengths we make use of the VLA 1.4-GHz 
source catalogue from Morrison et al. (2010) and Owen et 
al. (2008) for the GOODS-N and Lockman Hole fields, re- 
spectively. In GOODS-N, the Morrison et al. (2010) imaging 
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have a sensitivity of 3.9 pjy beam -1 with the sensitivity de- 
clining to 8pJybeam -1 at a radius of 15 arcmin from the 
centre. In Lockman Hole North, the Owen et al. (2008) imag- 
ing has a sensitivity of 2.7 pjy beam -1 over the central 40 
arcmin x 40 arcmin region overlapping the mm and SPIRE 
observations. 


3 MULTIWAVELENGTH IDENTIFICATIONS 

FOR SUBMM /MM DETECTED SOURCES 

The large beam size (> 10 arcsec FWHM) of typical 
submm/mm imaging facilities means that accurately com- 
piling multiwavelength identifications and photometry for 
submm/mm detected sources is often a significant challenge 
(e.g. Lilly et al. 1999; Ivison et al. 2007; Roseboom et al. 
2009; Chapin et al. 2011). This challenge arises mainly in 
two ways; firstly the large beam size, and low signal-to- 
noise, of typical submm/mm surveys means that catalogued 
sources will have quite large positional uncertainties. This 
makes matching to catalogues at other wavelengths difficult, 
as there may be more than one potential counterpart within 
the positional errors of the submm/mm source. Indeed, for 
faint submm sources at high redshift, the source density of 
matching catalogues deep enough to contain the true match 
may be so high that more than one match can always be 
found. 

The second problem is source confusion; more than one 
astronomical source is present in the submm/mm beam, and 
hence catalogued sources are made up of multiple astro- 
nomical objects. This problem is particularly worrying as 
it affects the integrity of the submm/mm source catalogues 
themselves. 

Several statistical methods which can account for one 
or both of these issues exist (Downes et al. 1986; Chapin et 
al. 2011; Roseboom et al. 2009; Roseboom et al. 2010). For 
the datasets considered here one or both of these issues are 
present. Source identification is a major issue for the mm- 
wavelength datasets, but source confusion is not, as mm- 
sources have large positional uncertainties, but low source 
densities. The SPIRE datasets have both large positional 
uncertainties and high source densities, and so are equally 
affected by both. 

Thus different approaches are needed to reliably 
produce multi-wavelength associations for the different 
submm/mm datasets. In practice this means that it is gen- 
erally not possible to reliably identify SPIRE counterparts 
to mm-detected sources without first making identifications 
in some other catalogue with higher spatial resolution (i.e. 
1.4 GHz or 24 /jm). To emphasise this point; the typical 
flux density of our mm-detected sample, after correcting for 
flux boosting, is Si. 2 mm ~ 2 mjy. Taking rough estimates 
of the mean dust temperature and redshift of submm/mm 
sources; 35 K (Kovacs et al. 2006; Chapman et al. 2010) and 
2 = 2.5 (Wardlow et al. 2011), and assuming a single tem- 
perature modified blackbody SED, gives predicted SPIRE 
fluxes of 13, 19 and 15 mjy at 250, 350 and 500 pm, re- 
spectively. These flux levels are below the confusion ‘limit’ 
at these wavelengths, with the best estimate of the noise 
from confused sources <r con f = 6mJy (Nguyen et al. 2010; 
Glenn et al. 2010). Hence recovering the SPIRE counter- 
parts to mm-detected sources will require a technique which 
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can reduce the effect of confusion, by using prior informa- 
tion at other wavelengths (Roseboom et al. 2010, Chapin et 
al. 2011). Putting the issue of source confusion aside, even 
if these sources could be recovered reliably, source identifi- 
cation would still pose a serious challenge. The positional 
uncertainties of both the mm and SPIRE detected sources 
(opos ~ 2-3 arcsec) requires the use of a 9-12 arcsec match- 
ing radius in order to recover the bulk of the true mm-SPIRE 
associations. Given known number density of SPIRE sources 
at these flux levels (~ 7000 deg -2 ; Glenn et al. 2010) we 
would expect between 14-25 per cent of these associations 
to be chance alignments. 

In this work we first match the mm-detected source 
lists to 1.4 GHz sources, using the well established catalogue 
based cross- matching techniques of Downes et al. (1986) 
and then use these 1.4 GHz source positions, along with the 
source positions of all known 1.4 GHz and 24 /rm sources in 
these fields, to ‘deblend’ the SPIRE maps into their indi- 
vidual contributors and thus provide the best estimate of 
the SPIRE photometry at the locations of our mm-selected 
sample. 

3.1 1.4 GHz identifications of mm-selected sources 

In the GOODS-N field we search for potential 1.4 GHz coun- 
terparts in the Morrison et al. (2010) VLA source catalogue 
to the 41 1.16-mm sources from Penner et al. (2011) us- 
ing a 10 arcsec search radius. 28 of the 41 sources have at 
least one potential counterpart within this search radius. To 
determine the reliability of these matches, we estimate the 
probability of a chance alignment using the P-statistic, as 
defined by Downes et al. (1986). Defining a reliable match 
to be one with a less than 5 per cent probability of being 
a chance alignment, 32 reliable matches are found for 24 
objects, with 8 sources having more than one reliable radio 
counterpart. 

For the sources with more than one reliable counterpart 
we take the counterpart with the lower P-statistic to be the 
correct identification. 

One difficult case for our matching is HDF 850.1 (ID 14 
in the Penner et al. 2011 catalogue), the brightest submm 
source detected in the original 850 gm SCUBA imaging of 
GOODS-N (Hughes et al. 1998). Our naive cross matching 
approach would associate this object with a 56 /jJy radio 
source located 5.3 arcsec away. However subsequent near-IR 
and high resolution Submillimctre Array (SMA) observa- 
tions show that the true position is much closer to the origi- 
nal SCUBA position (Dunlop et al. 2004; Cowie et al. 2009), 
although no distinct 1.4 GHz radio source can be found near 
this position in the Morrison et al. (2010) image. Here we 
adopt the Dunlop et al. (2004) position for matching with 
the SPIRE data. 

Other difficult cases are those sources for which the 
submm emission is later shown, via high resolution mm- 
wavelength interferometry, to originate from multiple com- 
ponents. In these cases the 1.4 GHz identification(s) is often 
found to lie co-incident with one (or more) of these com- 
ponents (e.g. Wang et al. 2011). Amongst our sample are 
two examples of this, GN20 (G1 in our sample; Daddi et 
al. 2009a; Carilli et al. 2010; Carilli et al. 2011) and GN21 
(G19 in our sample; Wang et al. 2011). In both cases we have 
simply taken a single radio identification, the one with the 


lowest P-statistic, as being solely responsible for the mm- 
flux density. For G1 the radio identification corresponds to 
GN20 from Daddi et al. 2009a, while for G19 this corre- 
sponds to GOODS 850-13c in the list of Wang et al. 2011. 
While it may turn out that a large fraction of mm-detected 
sources do indeed contain multiple components, in both G1 
and G19 the single radio identification we have chosen is 
found to be responsible for ^ 50 of the submm/mm flux 
(Daddi et al. 2009a; Wang et al. 2011). Importantly, our 
SPIRE photometry pipeline accounts for all 1.4 GHz, mm, 
and 24 /j,m detected sources, so it should not be affected by 
this issue. 

Encouragingly, for the 28 sources from Penner et al. 
(2010) which were previously identified in the AzTEC 
and/or MAMBO surveys we recover the same identifications 
made by both Chapin et al. (2009), for the 25 sources in com- 
mon, and Greve et al. (2008), for the 8 sources in common. 

In the Lockmann Hole North field we make use of the 
existing matches between the 1.2-mm sources and 1.4 GHz 
radio catalogue presented in Lindner et al. (2011). 40 of the 
41 1.2-mm sources are found to have a 1.4 GHz counterpart 
within 8 arcsec, with 39 of these found to have less than a 5 
per cent probability of being a chance alignment and hence 
are deemed reliable. 

It is worth noting that 3 of these reliable radio coun- 
terparts are not found in the published Owen et al. (2010) 
1.4 GHz source lists, and come from a deeper extraction per- 
formed in the vicinity of the MAMBO sources by Lindner 
et al. (2011). For these sources we adopt the quoted 1.4 GHz 
positions and flux densities from Lindner et al. (2011). 

Thus we are left with a sample of 63 mm-selected 
sources with high angular resolution 1.4 GHz positions with 
which to match to the SPIRE data. 

Interestingly there is a significant disparity between the 
identification rates in the 2 fields considered here. Both Pen- 
ner et al. (2010) and Lindner et al. (2011) list 41 mm- 
detected sources with a similar noise; la ~ 0.7mJy, al- 
though in both cases the noise varies significantly across 
the field. However the effective area of the two surveys 
is quite different; The Penner et al. (2010) combined 
AzTEC/MAMBO image of GOODS-N covers an area of 
0.08 deg. 2 , while the Lindner et al. (2011) MAMBO image 
of Lockman North covers an area of 0.16 deg. 2 . Thus it is 
clear that for the Penner et al. GOODS-N sample to find 
the same number of objects in similar to one half the area it 
must be identifying a fainter population of sources. The dif- 
ferential number counts at these wavelengths are known to 
be well described by a power law dN/dS oc S s , with an ex- 
ponent S ~ —3. So to observe a difference of a factor of 2 in 
surface density, the GOODS-N sample must be effectively a 
factor of \pl deeper, ignoring the fc-correction between 1.16 
and 1.2 mm. 

3.2 SPIRE Photometric method 

The most difficult obstacle to determining accurate SPIRE 
photometry is the effect of source confusion, i.e. the con- 
tributions of numerous faint sources within a single SPIRE 
resolution element, centred on the target of interest. For 
this reason, several authors have developed techniques that 
utilise the positions of sources detected at other wavelengths, 
usually 24 /.im and 1.4 GHz, to disentangle the various con- 
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tributions from discrete sources to the SPIRE flux in a 
given beam element (e.g. Chapin et al. 2010, Roseboom et 
al. 2010). This process is made possible by the high corre- 
spondence between the 24-pm and 1.4-GHz populations and 
those observed at Far-IR wavelengths; > 80 per cent of the 
cosmic IR background (CIB) at SPIRE wavelengths can be 
accounted for by 24-pm sources with £24 > 25 pjy (Pascale 
et al. 2009; Marsden et al. 2009), while the strong correlation 
between the Far-IR and radio luminosity is known to hold 
across a wide range in redshift and luminosity (e.g. Ibar et 
al. 2008; Ivison et al. 2010a). 

Here we follow the prescription presented by Roseboom 
et al. (2010; hereafter R10) to measure the SPIRE fluxes of 
mm sources in GOODS-N and Lockman Hole North, with 
a number of modifications. Below, we summarise the R10 
methodology, then highlight our departures from it. 

R10 assume that the SPIRE map can be Et directly 
by assuming that the positions of all sources contributing 
significantly to the map are known (i.e. previously detected 
at 24 pm, 1.4 GHz, or other wavelengths), and that only the 
SPIRE flux density of each of these sources is unknown. 
Algebraically this means the observed SPIRE map, d, can be 
described as a series of point sources with fluxes, /, at known 
positions, x. Thus, the observed map can be described as 

d = P x f + 8, 

Where P x is the point response function (PRF) centred 
at position x and 8 is some unknown noise term. As dis- 
cussed in R10, this equation can be easily inverted, using 
linear algebra, to obtain the best-fit values for the fluxes. 
However, one potential pitfall of this approach is that of 
over-fitting. If many of our input sources are intrinsically 
faint and the source density is high, this approach will give 
poor results. R10 solve this issue by using model-selection 
techniques to eliminate those sources that are unnecessary 
to describe the map. Specifically the Akiake Information 
Criterion (AIC; Akiake 1974) is used to filter the input 
list of sources via backward elimination. This is very time- 
consuming and hence the number of sources that could be 
fitted simultaneously was limited to ~100-200 in R10. 

Here, we replace this model-selection approach with the 
adaptive Least Absolute Shrinkage and Selection Operator 
(LASSO Tibshirani 1996; Zou 2006). In particular, we 
combine a non-negative application of this technique - the 
Non-Negative Least Squares path (NNLSpath) algorithm of 
ter Braak et al. (2010) - with the adaptive weighting scheme 
suggested by Zou (2006). In Zou (2006), the weights are 
tunable, based on initial estimates of the parameters; here 
we use the probability of a chance alignment, </, as defined 
in R10. Specifically, (j) is calculated separately for the three 
samples which make up our input list; radio, 24 pm and mm- 
detected sources. In cases where a source is detected in more 
than one of these bands the lowest value of <j> is adopted. As 
mm-detected sources have the lowest areal density these typ- 
ically have higher weightings, followed by radio sources, and 
then 24 pm. 

Accuracy in model selection is the primary advantage, 
the adaptive LASSO is known to have ‘oracle’ propertied] 
with the correct choice of weightings, however there is also 

3 In model selection, a method is said to have oracle properties 


a significant advantage in computational speed. With this 
approach we can solve for large numbers of candidate sources 
(up to 10,000) in a realistic amount of time. 

In addition, we modify the background estimation. 
Rather than use an unconstrained local estimate of the back- 
ground, the background is estimated globally, assuming a 
flat background across the total extent of the map. This ap- 
proach is prefered as it is more robust against incompleteness 
in the prior input list, and because there are no significant 
large scale correlations in the noise present in our SPIRE 
maps due to the exceptionally low frequency “knee” to the 
1 // noise (Griffin et al. 2010). 

Finally potential new SPIRE sources, i.e. those that 
were previously unknown at 24 pm or radio wavelengths, 
are considered by examining the residual maps after fitting 
the SPIRE fluxes of known sources. 

3.3 Estimating the true noise for SPIRE 
photometry method 

Accurately estimating the effect of confusion on an individ- 
ual source is non-trivial. As the SPIRE maps have a mean 
intensity of zero, for blind source catalogues we can expect 
that the typical confusion noise is given by the width of the 
distribution of flux densities (P(D)) in a PRF-convolved im- 
age. For SPIRE we know these values are 5.8, 6.3 and 6.8 
mjy/beam at 250, 350 and 500 pm, respectively (Nguyen 
et al. 2010). However our source extraction process both 
attempts to find the true zero point of the map, and simul- 
taneously fits the positions of sources found at other wave- 
lengths, in an attempt to reduce this noise. But it is difficult 
to assess for a given source how much we have ‘resolved’ the 
confusing background, and how much is left as unresolved 
fluctuations in the map. 

One simple way to assess the global improvement found 
by our method is to use completely synthetic maps and 
catalogues produced using input mock catalogues known to 
give reasonable agreement to the known far-IR and submm 
number counts (e.g. Fernandez-Conde et al. 2008). Testing 
against simulated SPIRE maps with the same properties as 
our GOODS-N data (i.e. instrumental noise, number counts 
and pixel intensity distribution) we find a typical total noise, 
i.e. confusion and instrumental noise, of ototai = 2.8, 3.9, 
3.7 mjy at 250, 350 and 500 pm, respectively. These values 
compare quite favourably to the SPIRE photometry process 
described in R10, which gave <7 to ta.i = 4.9, 5.4, 6.5 mjy at 
250, 350 and 500 pm for the same simulations. 

While testing the performance of our method on sim- 
ulated datasets is of interest, we really need an accurate 
estimate of the true error on the SPIRE photometry of 
our actual sources. One approach would be to estimate the 
global confusion noise in our SPIRE photometry by inves- 
tigating residual intensity fluctuations in the SPIRE maps, 
after subtracting the contributions from ‘resolved’ sources. 
This approach is appealing, as we know that the residual 
maps contain information about sources not included in our 
prior source list, as well as any artefacts from incorrect as- 
trometry and/or assumptions about the PRF. 

if it can reliably identify which free parameters are truly present 
in the underlying system, and which are unnecessary or spurious. 
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However there is one complication for this approach; 
variations in sensitivity in the prior source lists across 
the SPIRE map will mean the ability to resolve confusing 
sources varies spatially. Thus we must assess the confusion 
noise on scales small enough that the prior source list sensi- 
tivity is invariant. This is particularly difficult in this work, 
as the combination of 24 /jrn and 1.4 GHz prior source cata- 
logues results in large variations across the held. This is due 
to both the difference in relative sensitivity of 24 /jrn and 
1.4 GHz to SPIRE sources, as well as the tapering of the 
1.4 GHz sensitivity towards the edges of its coverage. 

We utilise a variant of this approach to estimate the 
true error for our confused SPIRE sources here. Specifically 
we construct PRF convolved residual maps for each of the 
SPIRE bands, using our best estimate of the SPIRE pho- 
tometry at the point of each input source and the known 
SPIRE beam. Then at the position of each source we are 
interested in we assess the total noise per pixel by measur- 
ing the standard deviation of surrounding pixels within a 20 
pixel radius (2, 2.8 and 4 arcmin for the SPIRE 250, 350 
and 500 /im maps used here). 

A further complication arises in converting this noise 
per pixel into the noise for a point source as we know that 
the confusing background must be correlated, on small scales 
by the PRF, and on larger scales by the intrinsic clustering 
of the sources. This means that in calculating the total noise 
estimate for a given point source we must take into account 
the uncorrelated instrumental noise in each pixel, and the 
correlated noise from confusion. For simplicity we assume 
that the correlated noise on small scales is dominated by 
the PRF, and assume the shape of the PRF gives the co- 
variance between pixels directly, i.e. If we assume the SPIRE 
photometry is estimated by 

E,W 

where di is the intensity in pixel i, Pi the PRF for the 
point source in pixel i and cn the instrumental noise in pixel 
i, the noise estimate, including contributions from both con- 
fusion noise and instrumental noise can be given by 

ffiotal = V (P T C- 1 P) , 

where C is the covariance matrix, the diagonal elements of 
which are given by 1/of, and the off diagonal elements i,j 
given by Pj/ai/aj, where P’ is the PRF assuming a point 
source centred at position i. 

Using this approach we find the typical ler total noise 
for our mm-detected sources is ototai =2.2, 2.6 and 2.5 mjy 
in the central part of the GOODS-N field for the 250, 350 
and 500 |im bands, respectively, increasing to 3.7, 4.8 and 
5.7 mjy at a radius of 20 arcmin. Similarly in the Lockman 
North field the total noise is found to be 3.5, 3.6 and 4.7 mjy 
at the centre of the field, and 4.4, 4.7, 6 mjy at a radius of 
20 arcmin. 

It is worth noting that these variations are solely due 
to the variations in the depth of the prior source list, the 
SPIRE sensitivity is uniform across both fields within the 
coverage of the mm observations (Levenson et al. 2010). 


4 RESULTS 

4.1 SPIRE detection statistics 

We define the SPIRE measurement as a detection if it has 
>3-cr significance in the flux density estimate and the good- 
ness of the photometric fit to the map is high (reduced 
X 2 < 5, in a 15 x 15 pixel region centred on the source 
position). The X' 2 cut is necessary to remove cases where 
bright or unaccounted for sources, or residual non-Gaussian 
noise (e.g. Galactic cirrus, or uncorrected thermal drifts) in- 
troduce signal that cannot be accounted for by our input 
source model. 

Table [Tj presents our SPIRE photometry for all mm- 
selected sample in GOODS-N, while Table[2]presents SPIRE 
photometry for mm-sources in the Lockman Hole North 
field. 

In GOODS-N, 17 of the 24 1.16-mm sources with good 
identifications are detected in at least 1 SPIRE band, while 
in Lockman North 30 of the 39 sources with identifications 
are detected in at least 1 SPIRE band. These figures give a 
marginally lower detection rate in GOODS-N (71 per cent) 
as compared to Lockman North (77 per cent), although this 
is consistent with the slightly deeper nature of the GOODS- 
N mm-imaging (See 1 )3.11) . 

As GOODS-N and Lockman Hole North are well- 
established surveys fields a large number of spectroscopic, 
and good quality photometric, redshifts are available. In 
GOODS-N nine sources have spectroscopic redshifts, either 
from targeted observations of SMGs and radio sources by 
Chapman et al. (2005, and in prep) or from the many mag- 
nitude limited redshift surveys of the field which have been 
conducted (see Barger et al. 2008 for a recent compilation of 
spectroscopic redshifts in GOODS-N). In the Lockman Hole 
North spectroscopic redshifts are available for four sources, 
mostly from mid-IR spectroscopy of 24 /jm selected sources 
(Fiolet et al. 2010), or WIYN spectroscopy of radio sources 
(Owen & Morrison 2009). Photometric redshifts are avail- 
able for an additional nine mm-sources in Lockman Hole 
North from the catalogue of Strazzullo et al. (2010). Here 
we only consider photometric redshifts which have a quality 
flag of “AA”, i.e. the estimated photo -2 error is less than 
0.2(1 + 2 ), from Strazzullo et al. (2010). In total good qual- 
ity redshifts are available for 22 of our parent sample of 63 
mm-sources. These redshifts are quoted in Tables [Q and [5] 
All of the sources with known spectroscopic and/or reli- 
able optical/near-IR photometric redshifts are detected by 
SPIRE. This is not particularly surprising, as those sources 
with redshift information from the optical/near-IR. wave- 
length range tend to be either more intrinsically luminous 
and/or at slightly lower redshift (Wardlow et al. 2011). 

To estimate the reliability of our detections, we test 
random positions in the SPIRE maps. Test positions are 
required to be >3arcsec from other sources in our input 
list, and within the region covered by our 24-/rm and 1.4- 
GHz catalogue. This is the typical minimum separation of 
sources in our input list, due to the resolutions at 24 /im 
and 1.4 GHz. If we impose the same selection criteria on the 
fluxes recovered for these random positions, 2.7 per cent are 
detected in any SPIRE band at a significance of >3<r 

Taking this as a guide, if our parent sample of 63 mm- 
sources were completely spurious then we would expect to 
recover SPIRE detections in any band for 2 sources. 
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Table 1. Measured SPIRE photometry for mm-sources in GOODS-N. 3a upper limits are given for sources which are not detected above 
this significance in the SPIRE imaging. For the mm photometry deboosted flux densities are quoted. The x 2 represents the reduced x 2 
in a 15 x 15 pixel region centred on the source position. ID number refers to the equivalent value from Penner et al. (2010). C09 refers 
to the ID numbers of sources in common with Chapin et al. (2009). G08 refers to the ID numbers of sources in common with Greve et 
al. (2008). 1.4 GHz flux density is from the catalogue of Morrison et al. (2010). 


Name 

C09+ 

G08* 

R.A. 

(deg) 

Dec. 

(deg) 

S250 

(mJy) 

v 2 

*250 

S 35 O 

(mJy) 

*350 

S 500 

(mJy) 

v 2 

*500 

.2mm 

(mJy) 

Si. 4GHz 
(mJy) 


Gl* 

1 

1 

189.29950 

62.37013 

18.6 ± 2.3 

0.97 

31.7 ± 2.7 

2.16 

24.3 ± 2.7 

1.10 

10.5=1=0.7 

0.079 

=t 

0.014 

G2* 

3 

2 

189.13918 

62.23579 

18.1 ± 2.3 

1.77 

15.2 ± 2.7 

2.11 

< 7.8 

1.68 

5. 2=1=0. 6 

0.036 

=t 

0.004 

G3* 

5 

4 

189.37811 

62.21633 

46.3 ± 2.0 

1.25 

46.5 ± 2.3 

1.75 

32.7 d= 2.3 

1.27 

4. 4=1=0. 6 

0.123 

± 

0.004 

G4* 

7 

3 

189.29727 

62.22540 

38.6 ± 2.1 

1.76 

33.0 ± 2.4 

2.62 

23.2 ± 2.3 

1.30 

4. 0=1=0. 5 

0.131 

± 

0.009 

G5 

2 

13 

189.13301 

62.28749 

< 6.9 

1.58 

< 8.7 

1.68 

< 7.8 

0.88 

4.1±0.6 

0.026 

=t 

0.005 

G6* 

6 

5 

189.11346 

62.10158 

9.2 ± 2.4 

1.23 

13.5 ± 3.0 

1.61 

5.2 ± 2.9 

1.14 

4. 5=1=0. 7 

0.053 

± 

0.012 

G7 

4 

12 

188.95979 

62.17839 

< 7.5 

1.53 

< 8.1 

2.31 

< 8.4 

2.09 

4.6±0.7 

0.028 

± 

0.005 

G8* 

26 

6 

189.30772 

62.30739 

8.2 ± 2.2 

1.91 

8.1 ± 2.8 

3.76 

< 7.8 

2.94 

3.1±0.6 

0.652 

± 

0.004 

G9* 

11 

14 

189.14937 

62.11888 

33.9 ± 2.1 

1.45 

20.3 ± 2.7 

1.16 

< 8.4 

1.43 

2. 9=1=0. 6 

0.046 

=t 

0.008 

Gil* 

13 

15 

188.97157 

62.22708 

21.0 ± 2.4 

1.25 

< 9.3 

4.06 

< 8.4 

2.39 

2. 6=1=0. 6 

0.045 

± 

0.005 

G13 

12 


189.13604 

62.10586 

< 6.9 

1.44 

< 8.7 

1.35 

< 8.7 

1.25 

2. 7=1=0. 7 

0.027 

=t 

0.005 

G14 

14 

25 

189.21696 

62.20714 

< 6.3 

1.05 

< 7.5 

1.77 

< 7.2 

1.31 

2. 4=1=0. 6 

0.016 

± 

0.004 1 

G16 



189.20117 

62.35201 

< 9 

1.00 

< 9.3 

1.98 

< 8.1 

1.24 

2. 5=1=0. 6 

0.030 

± 

0.005 

G18* 

16 


189.06708 

62.25388 

< 7.2 

1.35 

< 9.3 

2.95 

14.1 ± 2.6 

1.46 

2 .1±0.6 

0.036 

± 

0.005 

G19* 


29 

189.30020 

62.20341 

11.2 ± 2.1 

1.50 

11.3 ± 2.4 

1.28 

17.3 ± 2.4 

1.09 

2. 1=1=0. 6 

0.032 

± 

0.004 

G21* 

9 


189.40941 

62.29355 

7.1 ± 2.0 

1.10 

12.1 ± 2.6 

1.53 

< 8.4 

1.01 

2 .1±0.6 

0.026 

± 

0.006 

G24 

24 


189.03583 

62.24325 

< 7.5 

1.34 

< 9.6 

3.08 

< 8.1 

2.36 

2. 0=1=0. 6 

0.039 

=t 

0.004 

G26* 

15 


188.94971 

62.25811 

< 8.1 

2.33 

< 10.2 

4.22 

10.1 d= 2.9 

1.93 

2. 1=1=0. 7 

0.124 

± 

0.006 

G27* 

18 


189.42151 

62.20591 

31.0 ± 2.3 

1.19 

27.9 ± 2.4 

0.99 

30.3 ± 2.5 

0.88 

1. 9=1=0. 6 

0.033 

± 

0.005 

G29* 

25 


189.21550 

62.08417 

12.6 ± 2.9 

1.78 

15.1 ± 3.3 

2.57 

17.7 ± 3.0 

1.43 

2. 6=1=0. 9 

0.079 

± 

0.006 

G30* 

17 


188.92167 

62.24403 

11.7 d= 2.6 

1.59 

< 10.2 

3.80 

< 9. 

2.27 

2 . 2=1=0. 7 

0.029 

± 

0.005 

G31* 



189.32254 

62.13453 

35.1 ± 2.5 

3.20 

32.1 ± 3.0 

1.95 

22.5 ± 2.7 

0.96 

2 . 0=1=0. 7 

0.050 

± 

0.008 

G34 



189.14537 

62.32322 

< 7.2 

1.36 

< 8.7 

2.28 

< 10.3 

0.99 

1. 8=1=0. 6 

0.085 

± 

0.005 

G41* 



189.55217 

62.24859 

25.3 ± 3.2 

2.48 

25.6 ± 3.5 

9.72 

15.3 ± 3.7 

2.97 

1. 9=1=0. 7 

0.056 

± 

0.006 


i Typically written with prefix AzGN 
^ Typically written with prefix GN1200. 

*Deboosted flux density as described in Penner et al. (2011) 
* detected in at least 1 SPIRE band 
1 1.4 GHz flux density from Dunlop et al. 2004 


4.2 Submm-mm colours 

For the use of submm/mm photometry, either alone or in 
combination with mid-IR or radio data, to be useful as a 
redshift estimator there must be a reasonably strong, and 
unique, matching between submm/mm colour evolution and 
redshift. 

In Fig. □ we compare the SPIRE-mm colours to the 
expected evolution for a modified blackbody parameterised 
by a single dust temperature, Td, and optical depth r = 
(u/uo) heta , i.e. 

S„ oc [1 - exp(— (— ) p ]B v (Td), 
vq 

where B „ is the Planck function for a single dust tem- 
perature Td. In order to give a better representation to the 
mid-IR SED, we replace the blackbody SED on the Wien 
side of the SED with a power law, S u oc v a , where a is 
fixed to -2 (Blain et al. 2003; Conley et al. 2011). To ensure 
a smooth SED the power law is used at frequencies where 
d\ogS / d\ogv ^ a. As we have a limited number of photo- 
metric bands to constrain this model, we fix the dust opacity 
and optical depth to values of /3 = +1.8 and uq = 100 /im, 
respectively. We adopt this definition of a modified black- 
body throughout this paper. 

The sources with known spectroscopic, or reliable 
optical/near-IR photometric redshifts from Tables [3] and [I] 
are coloured coded by redshift. 

While the errors on the SPIRE and mm photometry are 
large, the modified blackbody model can sufficiently explain 
the observed submm-mm colours of the majority of SPIRE- 


mm sources. This is qualitatively in agreement with recent 
results from Kovacs et al. 2011, although they warn that 
more sophisticated models (in particular one which takes 
account of the range of dust temperatures likely to occur 
within a star forming galaxy) is needed to match shorter 
wavelength data (i.e. less than rest frame 30 /im). A modified 
blackbody redshift and temperature are completely degen- 
erate at z < 4, thus changing the assumed dust temperature 
will change the predicted submm-mm colour as a function 
of redshift, but it will not change the position of the model 
track in colour-colour space. Thus if the dust temperature 
were known the redshift could be inferred directly, or vice- 
versa. 

It can be seen that the highest redshift objects are 
found to have the ‘reddest’ submm/mm colours, i.e. large 
S'mm/S'spiRE or S 500 / S 35 o flux ratios, supporting claims that 
identifying red, or ‘cold’, SPIRE/mm sources is an effective 
way of isolating very high redshift submm galaxies (e.g. Pope 
et al. 2010). 

4.3 IR Luminosity to Dust Temperature Relation 
for SPIRE-mm sources 

As a significant fraction of SPIRE-mm sources have good 
redshift information (22 of 46 SPIRE detected mm-sources) 
we can investigate the IR luminosity to dust temperature 
relation directly for our sample of SPIRE-mm sources. 

For each source with a spectroscopic or optical/near-IR 
photometric redshift (excluding the very low -2 source L29) 
we fit a modified blackbody to the SPIRE and MAMBO 
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Table 2. Measured SPIRE photometry for mm-sources in Lockman Hole North. 3cr upper limits are given for sources which are not 
detected above this significance in the SPIRE imaging. For the mm photometry deboosted flux densities are quoted. The x 2 represents 
the reduced x 2 in a 15 x 15 pixel region centred on the source position. ID number refers to the equivalent value from Lindner et al. 
(2011). 1.4 GHz flux density is from the catalogue of Owen et al. (2008), as presented by Lindner et al. (2011). 


ID 

R.A. 

Dec. 

S250 


2 

*250 

*^350 

*350 

S500 


2 

*500 

S 1 2mm (deboosted) 

Si. 4GHz 



(deg) 

(deg) 

(mjy) 



(mjy) 


(mjy) 



(mjy) 

(mjy) 


LI* 

161.75083 

59.01878 

80.1 ± 

3.4 

0.92 

63.2 ± 3.8 

0.85 

40.6 ± 

4.4 

0.89 

0 c+0.6 
d,5 -0.6 

0.279 

± 

0.015 

L2* 

161.61192 

59.09578 

< 11.7 


1.07 

12.5 ± 4.0 

0.99 

< 13.8 


1.08 

0 O+0.7 
3 - 8 -0.7 

0.036 

± 

0.006 

L3 

161.63112 

58.84889 

< 10.8 


0.84 

< 11.1 

0.92 

< 13.2 


0.87 

q Q + 0'8 

3 - 8 -0.7 

0.459 

± 

0.025 

L4 

161.53050 

58.90389 

< 10.8 


0.83 

< 11.4 

1.03 

< 14.1 


0.79 

9 7+0.5 
Z --0.5 

0.030 

± 

0.004 

L5* 

161.85592 

59.06014 

32.8 ± 

3.9 

1.21 

44.5 ± 4.1 

1.08 

35.6 ± 

4.6 

0.92 

4.0+ 0 - 8 

-0.9 

0.051 

± 

0.006 

L6* 

161.66112 

58.93681 

52.4 ± 

3.6 

0.91 

51.5 ± 3.6 

0.75 

17.5 ± 

4.4 

2.17 

O q + 0.4 

0.160 

± 

0.010 

L7* 

161.75054 

58.91150 

21.0 ± 

3.6 

0.84 

18.6 ± 3.5 

0.97 

31.7 ± 

4.4 

0.78 

9 q+ 0.4 
Z '^ — 0.5 

0.042 

± 

0.005 

L8* 

161.63779 

58.86642 

30.0 ± 

3.6 

0.99 

29.8 ± 3.7 

0.94 

16.5 ± 

4.4 

1.25 

9 y+0.6 
— 0.6 

0.098 

± 

0.010 

L9 

161.77042 

58.83556 

< 10.8 


0.96 

< 11.1 

0.91 

< 13.8 


1.21 

3 - 8 ±i:S 

0.023 

± 

0.005 

L10* 

161.59604 

58.99300 

13.3 ± 

3.5 

0.72 

< 14.1 

0.94 

19.8 ± 

4.4 

0.91 


0.078 

± 

0.007 

Lll* 

161.48721 

58.88856 

58.2 ± 

3.7 

1.02 

56.6 ± 3.8 

1.26 

18.8 ± 

4.6 

1.44 

9 7+0.6 
-0.6 

0.315 

± 

0.019 

L12* 

161.19829 

59.00997 

25.3 ± 

3.7 

0.84 

36.0 ± 3.9 

1.11 

24.4 ± 

4.7 

0.78 

9 7+O.6 
- 0.7 

0.273 

± 

0.019 

L13* 

161.53646 

58.97458 

12.1 ± 

3.6 

0.98 

17.4 ± 3.6 

0.80 

< 14.1 


0.82 

2 l +0 - 5 

0.198 

± 

0.010 

L14* 

161.64937 

59.13014 

39.0 ± 

3.9 

0.99 

46.2 ± 4.0 

1.07 

33.9 ± 

4.6 

1.65 

3 q+0.9 
3 -0.9 

0.097 

± 

0.006 

L15* 

161.86654 

58.87058 

55.8 ± 

3.9 

1.01 

69.1 ± 3.6 

0.95 

49.7 ± 

4.8 

0.80 

3 n+0.9 
3-U -1.0 

0.180 

± 

0.018 

L16* 

161.83633 

58.86472 

< 11.7 


1.04 

11.2 ± 3.6 

0.90 

< 14.1 


0.95 

9 7+O.8 
-0.8 

0.048 

± 

0.007 

L17* 

161.54250 

59.04508 

10.4 ± 

3.7 

1.15 

14.6 ± 3.9 

1.26 

< 13.5 


0.84 

2 '°to .5 

0.027 

± 

0.003 

L18* 

161.73021 

58.83444 

39.9 ± 

3.6 

0.81 

37.8 ± 3.7 

0.91 

25.4 ± 

4.6 

0.95 

2-®±i:i 

0.105 

± 

0.008 

L19 

161.25833 

59.06758 

<11.4 


0.91 

< 12.3 

1.23 

< 14.7 


1.45 

2 o + 0 - 6 
-0.6 

0.029 

± 

0.004 

L21* 

161.37575 

59.11008 

< 11.7 


1.02 

12.0 ± 4. 

1.11 

16.0 ± 

4.9 

1.03 

1 g+ 0.5 
-0.5 

0.036 

± 

0.004 

L22* 

161.51763 

59.08039 

32.7 ± 

3.6 

1.06 

23.7 ± 3.9 

1.01 

< 14.1 


1.00 

2 0 +0 - 6 
^• u -o.6 

0.166 

± 

0.010 

L23 

161.67129 

58.89050 

< 11.1 


0.97 

< 10.8 

0.92 

< 13.2 


0.72 

1 7 + 0 - 5 

x - 0.5 

0.032 

± 

0.005 

L24* 

161.25279 

59.12600 

11.0 ± 

3.6 

0.92 

< 12.0 

0.76 

17.6 ± 

4.6 

0.79 

1 a+0-6 

1 ’ 8 — 0.5 

0.068 

± 

0.009 

L25 

161.41679 

59.06333 

< 12.0 


1.25 

< 12.6 

1.66 

< 15.0 


6.59 

1 7 +0 - 6 
1 -0.6 

0.024 

± 

0.007 

L27* 

161.76000 

58.85086 

17.4 ± 

3.5 

0.91 

19.5 ± 3.6 

0.71 

< 13.8 


1.16 

2 0 +0 - 8 
^• U - 0.9 

0.078 

± 

0.014 

L28 

161.58708 

58.90944 

< 10.5 


0.74 

< 11.1 

1.25 

< 13.5 


1.04 

c+ 0.5 

1 - 5 - 0.5 

0.025 

± 

0.006 

L29* 

161.48129 

59.15453 

140.5 ± 3.9 

2.47 

64.4 ± 3.9 

0.83 

29.1 ± 

4.8 

1.46 

1 8+ 0 - 8 
1 ’ 8 — 0.8 

0.307 

± 

0.043 

L30* 

161.29329 

59.06864 

< 11.7 


1.02 

13.3 ± 4.1 

1.28 

< 15.0 


1.59 

1 5+ 06 
i - t3 -0.6 

0.074 

± 

0.007 

L31* 

161.60379 

58.89644 

< 10.5 


1.04 

20.8 ± 3.7 

1.08 

13.1 ± 

4.3 

0.89 

1 5+ 0 - 6 
i-t> -0.5 

0.043 

± 

0.004 

L32* 

161.41579 

58.90692 

18.6 ± 

3.8 

1.07 

21.7 ± 3.9 

0.98 

17.6 ± 

5.2 

1.08 

1 7+O.8 
- 0.7 

0.047 

± 

0.004 

L33* 

161.39608 

58.84714 

28.3 ± 

3.7 

1.19 

32.5 ± 3.8 

0.96 

26.5 ± 

5.1 

1.13 

2 7+ 1 - 2 
- 1.2 

0.071 

± 

0.010 

L34* 

161.22346 

58.97819 

< 10.8 


0.94 

13.3 ± 3.8 

0.91 

< 14.1 


0.74 

1 5+ 0 - 7 
1 - 5 -0.7 

0.041 

± 

0.008 

L35* 

161.82546 

58.92383 

19.7 ± 

3.7 

0.95 

24.4 ± 3.7 

0.97 

18.2 ± 

4.5 

1.14 

1 5 +0 - 7 
i,b - 0.7 

0.060 

± 

0.006 

L36 

161.53375 

59.12889 

< 11.7 


1.40 

< 11.7 

0.98 

< 14.1 


1.13 

1.7+ 11 

-0.9 

0.016 

± 

0.003 

L37* 

161.54575 

58.87914 

52.0 ± 

3.6 

1.01 

39.9 ± 3.8 

1.09 

16.8 ± 

4.5 

0.95 

1 c+0.8 
1 - 5 -0.7 

0.160 

± 

0.011 

L38 

161.18704 

59.13836 

< 12.0 


0.86 

< 11.7 

0.81 

< 14.4 


4.57 

1 5+ 0-9 
i - t3 -0.9 

0.068 

± 

0.016 

L39* 

161.55012 

59.04261 

< 14.1 


1.12 

< 11.4 

1.22 

23.3 ± 

4.5 

0.81 

1 4+°- 7 

0.034 

± 

0.003 

L40* 

161.74304 

59.10931 

20.4 ± 

3.7 

1.09 

12.2 ± 4.0 

0.89 

< 14.1 


1.05 

1 5+ 1 2 

0.057 

± 

0.005 

L41* 

161.50129 

58.91794 

12.1 ± 

3.7 

0.86 

11.6 ± 3.8 

0.87 

< 14.1 


1.19 

1 4+°- 7 
1 ‘^-0.6 

0.072 

± 

0.012 


* detected in at least 1 SPIRE band 


photometry (and AzTEC and SCUBA where available), as- 
suming a dust emissivity of /3 = +1.8, the redshift and the 
functional form for a modified blackbody given in Section 
14.21 In the Lockman Hole North field PACS 100 and/or 
160 /im photometry is available from the HerMES survey 
for six sources and is included in the SED fits. In GOODS- 
N SCUBA 850 /jm photometry is available for five sources 
from the catalogue of Pope et al. (2005) and is included in 
the SED fits. IR luminosities are calculated for each source 
by integrating the modified blackbody fit in the range 8- 
1000 fj,m. IR luminosities and dust temperatures for sources 
with good redshift information in GOODS-N and Lockman 
Hole North are presented in Tables [3] and [4] respectively. 

Fig. [5] shows the best fit modified blackbody IR lumi- 
nosity vs. dust temperature for our sample of SPIRE-mm 
sources. One source, Gil, is omitted as our modified black- 


body model is not well-constrained by the observed data 
(See Table 0. Also shown is the local relationship for IRAS 
1.2 Jy sources from Chapman et al. (2003); after converting 
from IRAS Sioo/Seo colour to Td via our definition of a 
modified blackbody. It is clear that the SPIRE-mm sources 
form a tight sequence in Ljr-Td space. However we can 
see that our sample is systematically biased towards colder 
dust temperatures when compared to the local Ljr-Td re- 
lationship (Chapman et al. 2003). This bias toward colder 
dust temperatures for submm/mm selected samples is now 
well established (Magdis et al. 2010; Chapman et al. 2010; 
Magnelli et al. 2010). 

The bias towards colder dust temperatures causes the 
Lir Td relation to be much tighter than for an unbiased 
sample. To parametrise the observed Lir-Td relation for 
mm-selected sources we fit the functional form given by 
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Figure 1. Submm-mm colour evolution of mm-selected sources; 
*$ 500 / S 350 versus 535o/525o(t°p) ? 5'i.2mm/5'350 versus 5350/5250 
(middle) and 5i.2mm/55oo versus S 500 /S 350 (bottom). Only 
sources with > 3 <j detections in each of the relevant SPIRE bands 
are shown in each panel. The track taken by a modified blackbody, 
with Td = 40 K and optical depth r = (v/vq)@ where /3 = +1.8 
and vq = c/100 11 m is also shown. Sources located in the GOODS- 
N field are shown in squares, those in Lockman Hole North are 
shown as triangles. Symbols are colour-coded by redshift, with 
filled symbols designating those sources with spectroscopic red- 
shifts and open symbols those with reliable ( 5z < 0.2(1 + z)) 
optical/near-IR photometric redshifts. Black open symbols rep- 
resent sources with no redshift estimate. The submm-mm colours 
are reasonably well described by our simple modified blackbody 
model. 

Chapman et al. (2003) to our mm-selected sample. Specifi- 
cally Chapman et al. (2003) describe the mean IRAS colour 
via, 

log Co = log (7* + <51og(l + IfEL) +7 iog(i + (i) 

J^TIR J-^star 

where (7* = 0.45, S = —0.02, L* = 5 x 1O 1O L0, and 


Table 3. Radio-far IR photometric redshifts, far IR luminosities 
and dust temperatures for mm-detected sources in the GOODS-N 
field. Spectroscopic redshifts from the literature are quoted where 
available. Far IR luminosities (8-1000 fim) and dust temperatures 
are presented for sources with good redshift information. 


GOODS-N 


ID 

Redshift 

Radio/Far IR photo - z 

Lir(L©) 

Td(K) 

Gl* 

4.05“ 

+5.00+°;°° 

13.2±0.1 

39±1 

G2* 

4.05 6 

t 4 20 +0 ' 08 
^• zu - 0.21 

13.1±0.2 

43±2 

G3* 

2.98 c 

t 2 63" 1 " 0 ' 04 
^.UO— 0 06 

13.2±0.1 

50±2 

G4* 

1.996 d 

f 9 1 o+O.OO 

Z - icS -0.17 

12.7±0.2 

36±2 

G5 


*4.5 



G 6 * 


t 4 Q7+0-32 
-0.36 



G7 


*4.6 



G 8 * 


Aq 7+0.4 
' — 0.2 



G9 

0.95 e 

to 91+ 013 

u,yi — 0.15 



Gil* 

2.098 c 

1 0 QO + 0 - 22 
U,yZ -0.41 

12.8±1.3 

73±36 

G13 


*3.95 



G14 


*4.1±0.5 



G16 


*3.7 



G18* 

2.58 d 

t 2 50+O- 18 
Z. OO _0 27 

12.5±0.6 

35±8 

G19* 

2.91 d 

t2 91+ 0 - 24 
-0.23 

12.6±0.3 

36±3 

G21* 

3.19 d 

f 2.82t° 0 ;r 7 

12.6±0.4 

38±4 

G24 


*3.1 



G26* 


1 0 24+ 0 - 05 
u,z ^_ 0.09 



G27* 


t 2 43 + 0 - 04 
Z -^-0.18 



G29* 


t O ftl +0.26 

^•ol_0.31 



G30* 


t2 41+ 0,42 
0.26 



G31* 


*0.25+°'“ 



G34 


* 2.2 



G41* 


tl §7+0.17 
i -°‘ - 0.21 




* detected in at least 1 SPIRE band 

' Redshift estimate from simultaneous fitting of Lir-Td and 
LiR-Li. 4 Gifz relations 

t Redshift estimate from mm-1.4GHz spectral slope 
A Redshift estimate from fitting of Lir-Td as radio flux 
boosted by AGN activity. 

* Radio-Far IR photo- 2 : estimate from Dunlop et al. (2004) 
a Spectroscopic redshift from Daddi et al. (2009a) 

b Spectroscopic redshift from Daddi et al. (2009b) 
c Spectroscopic redshift from Chapman et al., in prep. 
d Spectroscopic redshift from Chapman et al. (2005) 
e Spectroscopic redshift from Barger et al. (2008) 


7 = 0.16, while Ltir is the integrated IR luminosity from 
3-1100 /im and Co is the ratio of the 60-to-100 /im IRAS 
flux density. We fit Eqn.|T|to our observed L/r-Td relation 
of our SPIRE-mm sources, adopting a conversion between 
IRAS 60-to-100 [i m flux ratio and Td given by a modified 
blackbody with /3 = +1.8. We leave L* fixed at the Chap- 
man et al. (2003) value as we do not have enough low lu- 
minosity sources to constrain this parameter. Our best fit 
is given by C* = 0.3, 5 = —0.03 and 7 = 0.193. The track 
taken by our best fit Ljr-Td relation is shown in Fig. [2] 
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Table 4. Radio-far IR photometric redshifts, far IR luminosi- 
ties and dust temperatures for mm-detected sources in the Lock- 
man Hole North field. Spectroscopic, or reliable ( 5z < 0.2(1 + z)) 
optical/near-IR photometric, redshifts from the literature are 
quoted where available. Far IR luminosities (8-1000 ^m) and dust 
temperatures are presented for sources with good redshift infor- 
mation. 


Lockman Hole North 

ID 

Redshift 

Radio/Far IR photo- £ 

Lir(Lq) 

Td(K) 

LI* 

2.562“ 

t 2 06 +0 ' 09 
z - uo -0.18 

13.3+0.2 

54+3 

L2* 


t4 QQ + 0 - 27 
-0.28 



L3 


* 1.8 



L4 


*4.4 



L5* 

3+0. 42^ 

t3 Q 0 +O - 07 

O.OU—o 29 

13.1+0.2 

43+2 

L 6 * 

2.03“ 

tl e 1+0-06 
± -° 1 -0.17 

12.9+0.2 

43+2 

L7* 


t 2 74+0 - 11 
z - '^- 0.20 



L 8 * 


+ 2 07 +o 2 ° 
z,u -0.23 



L9 


*3.9 



L10* 


t 3 0V +0.28 
°' U 1 - 0.21 



LIT 

1.95“ 

tl 71 +O.I 8 
A ' 1 -*--0.13 

12.9+0.2 

41+3 

L12* 

2 . 2 ± 0 . 2 i> 

t 2 64 -1 " 0 ' 08 
z,o ^-0.23 

12.7+0.2 

33+2 

L13* 


*0.32^ 



L14* 

2.3+0. 5 b 

*2.811°;!! 

12.9+0.2 

37+2 

L15* 

2 . 8 + 0 . 2 b 

t3 i 0+0.09 
o - ±o -0.16 

13.3+0.2 

45+2 

L16* 


t3 24+0- 51 
-0.33 



L17* 


+ 2 65” 1 ” 0 ' 29 



L18* 


t 2 15+0-H 
z,lo -0.24 



L19 


*4.1 



L21* 


f 2-52 



L22* 

1.4+0. 4 b 

f 0.38 ±g;|| 

12.3+0.5 

32+4 

L23 


*3.6 



L24* 

3.2+0. 4 6 

t 2 52+ 0 - 43 
Z, ' J - 0.20 

12.7+0.5 

40+5 

L25 


*3.5 



L27* 


tf 62 +0 ‘ 45 
I.OZ— o 27 



L28 


*3.8 



L29* 

0.004 d 

t 0.38±g;i| 

8 . 0+0.3 

22+2 

L30* 


to -7llS'.4? 



L31* 

2.9+0. 3 b 

t2 44 + 0 - 21 
0.24 

12.7+0.4 

40+4 

L32* 

1.3+0. 2 b 

t O 01 +0.25 
— 0.27 

12.0+0.4 

24+2 

L33* 


t 2 63 +0 16 
z -°°-0.34 



L34* 


t2 46 +0,50 
z -^°_0.32 



L35* 


+ 2 14+0-1 8 
z - 1 ^-0.24 



L36 


*4.5 



L37* 

1.7±0.2 fc 

00 

6 6 
+ 1 
L- 



L38 


*3.6 



L39* 


*2.591!;!! 



L40* 


t 0.78lg;i| 



L41* 


tf 4Q+ 101 
J -- i± »_0.87 




* detected in at least 1 SPIRE band 

f Redshift estimate from simultaneous fitting of Lir-T d and 
Lir-Li. 4 GHz relations 

J Redshift estimate from mm-radio spectral slope as presented 
by Linder et al. (2011) 

“ Spectroscopic redshift from Polletta et al. (2006) 
b Photometric redshift from Strazzullo et al. (2010) 

“ Spectroscopic redshift from Fiolet et al. (2010) 
d Spectroscopic redshift from Owen & Morrison (2009) 



fOQidiR (*- ) 


Figure 2. IR, luminosity vs. dust temperature for those SPIRE- 
mm sources with good redshift information. Source Gil is omit- 
ted as a modified blackbody model is not well-constrained by 
the observed data. Symbols are colour coded by redshift, with 
GOODS-N sources denoted by squares and Lockman Hole North 
sources denoted by triangles. The solid line is the local Lir-Td 
relation from the 1.2 Jy IRAS sample as determined by Chapman 
et al. (2003), where we have converted from .S’ i oo /Sqo to Tr as- 
suming a modified blackbody with optical depth r = (u / uq)^ 
with 0 = +1.8 and //q = c/100 (im. The dashed line shows 
a modification to this relation given by fitting Equation [T] di- 
rectly to our SPIRE-mm sources. The best fit parameters are; 
7 = 0.193, 5 = —0.03 and C* = 0.3. L* is held fixed at 5 X 10 10 . 

5 DISCUSSION 

5.1 Submm photometric redshifts 

As shown in 1 )4. 21 there is a reasonable correspondence be- 
tween submm/mm colour and redshift, which is likely a by 
product of our SPIRE-mm sources demonstrating a strong 
relationship between dust temperature and luminosity, as 
seen in Fig. [2] Thus, if we can assume a relationship be- 
tween dust temperature and IR luminosity, we can break the 
temperature-redshift degeneracy and produce photometric 
redshifts which rely on submm/mm data alone. 

We explore the potential of such an approach here by 
making use of the 22 SPIRE-mm sources with good redshift 
information as a test set. We assume that the Lir-Td rela- 
tionship for SPIRE-mm sources is given by our best fit of 
Eqn.[l]as seen in Fig. [2] Additional photometry from PACS 
and SCUBA is not considered in the photometric redshift es- 
timation as these are only available for a small subset of the 
sample; we want a fair assessment of the ability of SPIRE 
and mm-wavelength data alone to estimate the redshift. 

Fig. [3] compares our photometric redshift estimates to 
the known redshifts for these sources. The mean photomet- 
ric redshift accuracy is found to be |As|/(l + a) = 0.16 with 
the mean offset (< |Az| >= 0.51). All of our SPIRE-mm 
sources are detected in the radio at 1.4 GHz so we can hope 
to improve on our photometric redshift estimates by also 
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Figure 3. Radio-to-Far IR correlation as a function of redshift 
as probed by the giR parameter, where ^ir = log 10 [(L/R/3.75 x 
10 12 W)/(L L4 GJ i z / W Hz - 1 )] . The dashed line and shaded region 
represent the median and =bl a values for SPIRE sources from Ivi- 
son et al. (2010). Squares and triangles represent sources with reli- 
able redshift information in GOODS-N and Lockman Hole North, 
respectively. Good agreement is seen with the Ivison et al. (2010) 
values. No significant evolution with redshift is observed, suggest- 
ing our single qm model should introduce no redshift dependant 
bias to the photometric redshift estimates. 

requiring that the sources lie on the radio-to-Far IR cor- 
relation (van der Kruit 1971). As our SPIRE photometry 
is close to the peak of the Far-IR emission at z ~ 1 — 3, 
it is not sufficient to assume a spectral index between the 
SPIRE/mm bands and 1.4 GHz as is typically done for 
submm-radio photo - z (i.e. Carilli &; Yun 1999). Here we as- 
sume the ratio of Lir to Liaghz is given by a constant value 
qm = log 10 [(Lifl/3.75 x 10 12 W)/(Li. 4 gW W Hz' 1 )] = 2.4 
(Ivison et al. 2010). We also assume that the radio SED is ad- 
equately described by a power law, S„ oc v a , with a = —0.75 
(Ibar et al. 2010). Fig. [3] shows q as a function of redshift 
for the 22 SPIRE-mm sources with good redshifts. Good 
agreement is seen between our sample and the median (and 
standard deviation) of q for SPIRE sources as measured by 
Ivison et al. (2010). Importantly no evolution with redshift 
is observed, similar to Ivison et al. (2010), and hence our 
single q model should introduce no redshift dependant bias. 
As none of the sources in Fig. QT|lie significantly off the radio- 
to-far IR correlation, we are confident that potential contri- 
butions to the radio flux density from AGN activity may be 
ignored. 

To get the best estimate of the redshift we simultane- 
ously fit the Lir-Td and Lir-Li.4g.hz relations, assuming 
the Far-IR emission is described by a modified blackbody. 
After including the radio-to-Far IR correlation the photo- 
metric redshift accuracy improves to |Az|/(l + z) = 0.15 
or (< \Az\ >= 0.45). If we consider only those sources 
with spectroscopic redshifts the accuracy improves again to 
|Az|/(l + z) = 0.12 or (< |Az| >= 0.35). Interestingly this 
same improvement is not seen for the Lir-Td only method 
when the sample is restricted to only those with spectro- 
scopic redshifts (|Az|/(l + z) = 0.16 or (< \Az\ >= 0.55)). 

While finding the redshift that gives the maximum like- 


lihood redshift (i.e. the best fit) is the primary aim, it is also 
of interest to examine the probability distribution to see how 
wide a range in redshift also gives a reasonable fit. For each 
source we calculate the redshift range in which 68 per cent of 
the probability is enclosed. These points are shown as error 
bars in Fig. [I] Using this approach, the typical estimated 
redshift error for photometric redshifts using only the Lir- 
Td is a z = 0.6, while including the radio-to-Far IR correla- 
tion reduces this to cr 2 = 0.4. Encouragingly these numbers 
are quite similar to the true error measured by comparing 
to known redshifts. 

Many previous works have attempted to use a com- 
bination of radio and and ground-based submm/mm data 
to estimate photometric redshifts. Aretxaga et al. (2007) 
estimated submm-radio photometric redshifts for 120 850- 
/im sources detected in the SHADES survey (Mortier et al. 
2005), 58 of which had previously known spectroscopic or 
robust optical/near-IR photometric redshifts. They make 
use of two distinct photo -2 techniques; a simple 850 /im to 
1.4 GHz spectral index estimator, similar to that proposed 
by Carilli & Yun (1999), and full fitting of the Far-IR to ra- 
dio SED to a set of well-defined templates. Comparing the 
results of these two techniques to the known redshifts gives 
a mean photometric redshift error of | Az| = 0.9 for the spec- 
tral index method, and |Az| = 0.6 for the template fitting 
method. 

More recently Biggs et al. (2011) used a similar spectral 
index approach to estimate photometric redshifts for 870- /im 
sources detected in the LESS survey of ECDFS. Comparing 
the 48 sources which have spectroscopic redshifts, or reliable 
17-band photometric redshifts (Wardlow et al. 2011), the 
photometric redshift accuracy is found to be | Az|/(1 + z) = 
0.25 or a mean deviation of < |Aa| >= 0.68. 

Thus it is clear that photometric redshifts estimated 
using the Lir-Td relation alone offer a slight improvement 
over radio-submm redshifts, while using a combination of the 
Lir Td relation, and radio-far IR correlation, the accuracy 
of submm/mm-radio photometric redshifts can be increased 
by almost a factor of 2. By contrast, a similar study by 
Dannerbauer et al. (2010) using PACS 100 and 160 (im, in 
combination with mid-IR, radio, and submm/mm (but not 
SPIRE), data find that the accuracy of photometric redshifts 
using IR templates is significantly worse if PACS data is 
included. This highlights both the need for SPIRE data, 
which is at or close to the peak wavelength of the far IR 
emission at these redshifts, as well as the need for improved 
SED models in the mid IR (5-30 /im.) 

5.2 Redshift distribution of submm/mm sources 

The redshift distribution of sources selected at different 
submm and mm wavelengths has been a subject of much 
debate. Many previous works have claimed that sources se- 
lected at mm wavelengths are preferentially found at higher 
redshifts than those at ~ 850 gm (Dannerbauer et al. 2004; 
Younger et al. 2007, Greve et al. 2008, Coppin et al. 2009). 
In contrast, others have claimed that the ~850 gm and mm- 
selected populations are at similar redshifts, with the appar- 
ent differences originating from selection effects (Wardlow et 
al. 2011). 

Thus it is of interest to re-examine this using our sam- 
ple of SPIRE-mm sources. While only 22 of the 46 SPIRE 
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Figure 4. Spectroscopic vs photometric redshift estimates for submm/mm photometry alone, assuming a universal Ljr-Td relation, 
(left) and combined with radio constraints, assuming a universal radio-to-far IR correlation (right). Filled symbols designate sources with 
spectroscopic redshifts, while open symbols designate those sources with optical/near-IR photometric redshifts only. For submm/mm 
photometric redshifts we achieve a mean photo - z accuracy of \A.z\/ (l-\- z) = 0.16 (< |Az| >= 0.51), while if radio constraints are included 
this is reduced to | Ajz|/ (1 + z) = 0.15 (< |Az| >= 0.45). Error bars are calculated by finding the redshift range which contains 68 per 
cent of the probability. 


detected mm-sources have reliable spectroscopic or photo- 
metric redshifts, as we have seen in it is possible to 
estimate the redshifts of SPIRE-mm sources with reason- 
able accuracy using a combination of submm/mm and radio 
data alone. 

As the submm photometric redshift errors are large, and 
our sample of SPIRE-mm sources small, we prefer to use the 
full probability distribution function (PDF) for each source 
rather than simply taking the maximum likelihood estimate 
of the photometric redshift fit in constructing the redshift 
distribution. This has the advantage of giving a more robust 
reconstruction of the redshift distribution (Sheth & Rossi 
2010), although it should be remembered that the distribu- 
tion produced in this way is essentially a convolution of the 
true redshift distribution and the typical photo -2 error. 

In addition to the SPIRE-mm sources we can also esti- 
mate the redshifts of the 16 SPIRE undetected mm-sources, 
using the well-established radio-mm spectral index (Carilli 
& Yun 1999). For the sources in Lockman Hole North Lind- 
ner et al. (2011) quote estimates of the redshift for all 39 
1.2-mm sources with radio identifications, using, where pos- 
sible, a combination of 1.2-mm, 1.4 GHz and lower frequency 
GMRT 325 and 610 MHz observations. We use their redshift 
estimates for the nine sources without SPIRE counterparts 
in the Lockman Hole North field. In GOODS-N we again 
use the spectral index approach as described by Carilli & 
Yun (1999) to estimate redshifts for the seven SPIRE un- 


detected mm-sources, assuming the same values for a ra dio 
and a su bmm as Lindner et al. (2011), a ra dio = —0.68 and 
ctsubmm = 3.2, to ensure consistency between the two sam- 
ples. The one exception to this is G14, otherwise known as 
HDF850.1, for which we used the radio-Far IR photomet- 
ric redshift of 2 = 4.1 ± 0.5 determined by Dunlop et al. 
2004. To assess how reliable the mm-radio spectral index 
technique is for this sample we compare redshift estimates 
produced via this method to the known redshifts for the 
22 sources in our sample with good redshift information. 
We find a photo -2 accuracy of | A 2 1/( 1 + z) — 0.24, or a 
mean difference < | A 2 1 >= 0.6, for photometric redshifts 
estimated from the radio-mm spectral index alone. 

The full set of SPIRE-mm-radio and mm-radio pho- 
tometric redshifts is presented in Tables 0 and [I] for the 
GOODS-N and Lockman Hole North fields, respectively. 

The redshift distribution for our SPIRE-mm sources 
and SPIRE undetected mm-sources is shown in Fig. [5] 
SPIRE-mm sources show a strong peak at 2 ~ 2-3, simi- 
lar to other known submm/mm samples. Interestingly, the 
redshift distribution for SPIRE undetected sources shows 
that they are predominately found at 2 > 3. Only four of 
the SPIRE detected sample are estimated to be at 2 > 3.5, 
while 13 of the 16 SPIRE undetected mm-sources are found 
above this redshift. This is not surprising as we have already 
found in H4.2l that the sources with the largest mm-to-SPIRE 
colours are found to be at the highest redshifts. 
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Figure 5. Left: The redshift SPIRE-mm sources as determined by use of srBSfll^rlifti-radio photometric redshifts. 

SPIRE-detected mm-sources are shown as a grey shaded histogram. The distribution of SPIRE-mm sources with spectroscopic redshifts 
or reliable optical/near-IR photometric redshifts from Tables [l]and[2]is shown by the dashed line. Redshift estimates from the radio- 
mm spectral index for sources from our parent mm-sample without SPIRE detections are shown as a solid line. Right: The redshift 
distribution of SPIRE-mm sources (grey shaded) and all mm-sources (grey hatched) compared to 850 pm selected sources (solid blue 
line) from a compilation by Chapman et al. (in prep.), 1.1-mm selected sources (dotted red line) from Chapin et al. (2009.) and 870 pm 
sources (dashed line) from Wardlow et al. (2011). In each case the histogram is scaled to the total number of sources in our sample. The 
mm-selected sources are seen to peak at higher redshifts than those selected at ~ 850 pm, with a significant excess at z > 3 made up of 
primarily mm-only (i.e. SPIRE undetected) sources. 


However this large population of SPIRE undetected 
mm-sources at high redshift is somewhat at odds with the 
redshift distributions of previous submm surveys, which are 
not limited by SPIRE. Shown for comparison is the red- 
shift distribution of 850 pm selected sources, from a com- 
pilation by Chapman et al. (in prep.), 1.1-mm sources in 
GOODS-N from Chapin et al. (2009; henceforth C09), and 
870 pm sources in the LABOCA Extended Chandra Deep 
Field South Submm Survey (LESS; Weifi et al. 2009) from 
Wardlow et al. (2011; henceforth Wll). Comparing the red- 
shift distribution of SPIRE-mm sources and the submm 
(850/870 pm) selected samples of Chapman et al. (2005) and 
Wardlow et al. (2011; Wll) we find a clear deficit of z < 2.5 
SPIRE-mm sources. In addition, once the mm-only sources 
are included we find many more 2 > 3.5 sources than ei- 
ther Chapman et al. or Wll. Much better agreement is seen 
with the C09 sample, although this work shows somewhat 
fewer sources at low-z (i.e. z < 1), potentially due to the 
smaller area probed in GOODS-N compared to our larger 
sample. A K-S test suggests that the C09 redshifts and our 
mm-selected sample have a > 99 per cent chance of being 
drawn from the same underlying redshift distribution. 

One potential reason for the disagreement between the 
redshift distribution of the mm and 850/870 pm-sources is 
the depth of radio data available; in this work we have used 
the two deepest blank field images of the sky at 1.4 GHz 
available, with a typical rms noise of ~ 3-4 pjy. This is crit- 
ical as of the 13 a > 3.5 SPIRE undetected mm-sources, nine 
are fainter than 30 pjy at 1.4 GHz. C09 also make use of the 


Morrison et al. (2010) radio data in GOODS-N. In ECDFS 
the 1.4 GHz imaging has an rms of 6.5 pjy at its deepest, 
and thus the Wll will miss the counterparts for many of the 
highest redshift sources. Similarly the Chapman et al. list is 
based on lower quality radio data over several fields, and 
thus is also unlikely to identify the highest redshift submm 
sources. 

Wll attempt to account for submm sources which do 
not have robust optical/near-IR or radio counterparts by 
investigating the redshift distribution of optical-near IR 
sources neighbouring their submm detected sample, and this 
correction is included in the redshift distribution shown in 
Fig- El However while this means that the Wll redshift dis- 
tribution is not limited by the depth of the radio (or mid-IR) 
imaging used for identification, this correction is limited by 
the depth of the available optical/near-IR imaging. Wll es- 
timate that even after this correction, 20 per cent of LESS 
sources are still unaccounted for. 

While it appears that our work here confirms that radio- 
detected, mm-selected samples appear at higher redshift 
than those selected at ~ 850 pm, it is difficult to determine 
whether this result is not simply a product of our use of 
deeper 1.4 GHz radio data for identifications. To investigate 
the possibility that the A'-correction between submm (i.e. 
850 pm) and mm-wavelengths (i.e. 1.2-mm) is the cause for 
the difference in the redshift distribution in Fig. [6] we plot 
the limiting IR luminosity assuming a single modified black- 
body with dust temperature 30K and 40K as a function of 
redshift. The flux limits at 350 pm, 850 pm and 1.2 mm are 
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taken to be 10, 5, and 2mJy, respectively. The 350 /jm and 

1.2-mm limits were chosen to be representative of the sam- 
ple under study here, while the 850 /jm limit was chosen to 
be representative of current generation surveys at, or near, 
this wavelength (Weifi et al. 2009; Mortier et al. 2006). 

From Fig. [S] it is clear that, irrespective of dust tem- 
perature, 850 /im is marginally more sensitive in the red- 
shift range 1 < a < 3, while 1.2-mm is more sensitive at 
2 > 3. We can also see that the 350 /im limiting luminos- 
ity increases rapidly as a function of redshift, crossing the 

1.2-mm limit at z ~ 3. This explains the lack of SPIRE-mm 
sources seen above this redshift. 

While the difference in luminosity is small, the shape 
of the luminosity function at these extreme luminosities is 
very steep (Chapman et al. 2005; Wll). At 1 < 2 < 2 the 
850 /im and 1.2-mm sensitivities differ by ~ 0.1 decades in 
Ur- Using the Chapman et al. (2005) determination of the 
IR luminosity function, extrapolating to the limiting lumi- 
nosity shown in the left panel of Fig. [6] for a Td =30K SED 
(L/ij = 12.4 L©; log 10 ~ — 5.3 Mpc 3 dec. -1 ) and assuming 
the full volume of our parent surveys between 1 < 2 < 3 
(0.24 deg. 2 ; 5.4 x 10 6 Mpc 3 ) we could miss ~ 5 850 |im de- 
tected sources via a 1.2-mm selection. Similarly in the red- 
shift range 3 < 2 < 5, assuming the sensitivities differ by 
~ 0.2 decades in L ir and again using the Chapman et al. 
(2005) luminosity function estimate at Lir = 12.4 Lq, we 
could detect an extra ~ 10 sources via a 1.2-mm selection 
at the quoted flux limits. These numbers are crudely in line 
with the excess, and decrement, of sources we observe in red- 
shift distribution of our mm-selected sources at 1 < 2 < 3 
and 2 > 3, respectively. Performing a similar exercise assum- 
ing Td =40K would give a decrement of ~ 1 and excess of 
~ 2 sources at 1 < 2 < 3 and 3 < 2 < 5, respectively. Thus 
if the evolution of the 850 /rm-1.2 mm K correction is the 
cause of the observed differences in the redshift distribution 
it must be relatively cold (Td ^30K) sources which make 
up the discrepent population. 

Of course these calculations assume a particular optical 
depth t = (v/vo)P , although it is worth noting that the 
differences in sensitivity between 850 /im and 1.2 mm are 
enhanced for lower values of /3, and are relatively insensitive 
to Vo. 

5.3 Contribution of mm-selected sources to the 
SSFRd 

Now that we have some estimate of the redshift for all of 
our parent of radio detected mm-selected sample, we can es- 
timate their contribution to the IR luminosity density (pm) 
and hence the star formation rate density of the Universe 
(SFRd). Given the wide redshift distribution of mm-sources 
we estimate pm in four redshift bins; 0.5 < 2 < 1.5, 1.5 < 
2 < 2.5, 2.5 < 2 < 3.5 and 3.5 < 2 < 5. The redshift 
bins were chosen so as to match previous measurements 
(i.e. C05, Wll), while still maintaining a reasonable num- 
ber of sources per bin. The contribution of each source is 
estimated by taking the measured IR luminosity (Lir) and 
dividing by the maximum volume a source could be detected 
in (Unax). Incompleteness in the mm-sample is taken into 
account by dividing each object by the completeness at its 
mm-wavelength flux density. Completeness estimates for the 

1.2-mm source list in Lockman North are given by Lindner et 


Table 5. Contribution of mm-selected sources to the IR luminos- 
ity density. 


Redshift range 

PIR 

PSFR 


Logio Lq Mpc -3 

Logio M@ yr — 1 Mpc -3 

0.5 < 2 < 1.5 

7.1L0.1 

-2.7L0.1 

1.5 < 2 < 2.5 

7.72L0.07 

-2.05i0.07 

2.5 < 2 < 3.5 

7.86L0.08 

-1.91i0.06 

3.5 < 2 < 5 

7.7T0.1 

-2.1i0.1 


al. (2011). Unfortunately Penner et al. (2011) do not provide 
completeness estimates for their source list in GOODS-N, so 
we calculate the completeness as a function of 1.16- mm flux 
density by injecting sources at a range of flux densities into 
their combined MAMBO-AzTEC image. Sources are con- 
sidered recovered if they are found within 3 arcsec of the 
injected position with a peak signal-to-noise of a > 3.8. 

For each source we wish to use the most reliable red- 
shift available. For the 22 sources with spectroscopic or good 
optical/near-IR photometric redshifts we use these measure- 
ments, while for the 25 SPIRE-mm without prior redshift 
estimates we use our submm/mm-radio photo- 2 ’s from i l5.ll 
For the remaining 16 mm-only sources we use estimates from 
the radio-mm spectral index, as described in 35.21 

However, given the large photometric redshift errors 
of the submm/mm-radio and radio-mm spectral index es- 
timates which form the bulk of our sample we estimate 
pm via a Monte-Carlo approach in which realisations are 
producing by placing each source at a redshift randomly 
drawn from the probability distribution function (PDF) of 
the submm/mm-radio fitting. In each of the realisations 
sources with good spectroscopic redshifts are held fixed. 
Sources with optical/near-IR photometric redshifts are as- 
sumed to have Gaussian PDFs with the width determined by 
the quoted errors from Table [2] Similarly mm-only sources 
with radio-mm spectral index redshifts are assumed to have 
Gaussian PDFs with width given by the typical A 2 found 
in 35.21 A 2 = 0.6. In addition mm-only sources are assumed 
to be at a dust temperature of 35K. 

In each Monte-Carlo realisation we assess the contribu- 
tion of each source to each redshift bin by taking the sum 
7/" -I/iR/Vj/ax, where sum n runs over the sources which 
happen to fall in that redshift bin in a given realisation, and 
Lir and U ma x for source i are calculated using the best fit 
modified blackbody. To calculate Unax, we find the mini- 
mum and maximum redshift at which an observed source 
would still be found in our sample, i.e. > 3.8a at 1.16 or 

1.2-mm, depending on the field, and > 5a at 1.4 GHz, tak- 
ing into account variations in both the luminosity distance 
and K correction. The posterior probability of piR and the 
mean redshift in each bin is then inferred from a histogram 
of the Monte-Carlo realisations, fitting a Gaussian to return 
the peak likelihood and variance. 

Based on this analysis, our best estimate of the contri- 
bution of mm-selected sources to the IR luminosity density 
is given in Table [5] 

Fig. □ compares these values to previous measures of 
the contribution to the IR luminosity density by 850 pm 
and LESS 870 pm sources (Chapman et al. 2005; Ward- 
low et al. 2011), and a prediction of the contribution from 
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Figure 6. Limiting IR luminosity as a function of redshift for flux limited selection at 3 submm/mm-wavelengths; S 350 > lOmJy (solid 
line), .S'sso > 5mjy (dotted line) and Si. 2 mm > 2mjy (dashed line). In each case we have assumed a modified blackbody SED with dust 
temperature To = 30 K (left) and To = 40 K (right). In both panels, 850 //m is more sensitive in the redshift range 1 < z < 3, while 
1.2- mm is more sensitive at 2 > 3. 


Lm > 10 12 ' 5 Lq from the model of Bethermin et al. (2010). 
At 2 ~ 2.5 our estimates appear lower than both the con- 
tribution from ~ 850 gm sources, and the predictions for 
Lir > 10 12S Lq sources. However there are a number of 
ways to account for this. Firstly, as we have seen in 115.21 mm- 
selected surveys are less sensitive than their 850 /im counter- 
parts at 1 < z < 3. However mm-surveys must be less sensi- 
tive in general as the areal density of mm-detected sources 
is much lower than equivalent 850 gm samples. Including 
sources without radio identifications, our parent sample con- 
tains 82 mm-detected sources in 0.24 deg. 2 , while the LESS 
survey detect 126 sources at 870 fim in a 0.25 deg. 2 and 
SHADES detected 120 sources at 850 /im across two fields 
totalling 0.22 deg. 2 . 

Given this, it is puzzling that our estimate at 2 ~ 4 
is somewhat larger than that found for 870 fim sources by 
Wll, however as we have seen in 95.21 there are good rea- 
sons to expect a significantly larger contribution at z > 3, 
as mm-surveys are more sensitive at these redshifts than 
their 850 / tm counterparts, and we have the advantage of 
significantly deeper radio data with which to identify our 
mm-sources. The Wll analysis is missing the 45 per cent 
of LESS SMGs that do not have identifications, similarly 
Chapman et al. (2005) only have radio identifications for 
~70 per cent. Here we have radio identifications for 78 per 
cent of our parent mm-detected sample. 

It is of interest that our measurement of the SFRd 
at 2 ~ 4, while higher than the Wll estimate for 870 pm 
sources, is still somewhat short of that predicted for IR lu- 
minous galaxies (Logio Lir > 12.5Lq) by Bethermin et al. 
(2010), and ~ 100 lower than the total SFRd estimated from 
a combination of UV, radio and IR measurements (Hop- 
kins & Beacom 2006). Thus while selecting sources at mm- 
wavelengths offers an effective way to isolate luminous star- 
bursts at high redshift (z > 3), it is clear that mm-selected 


sources, at least at the imaging depth possible with current 
facilities, do not make a significant contribution to the SFRd 
at any redshift. 

Due to the strong dust temperature bias, mm-selected 
sources never contain the bulk of IR luminous sources at any 
redshift. Taking the Bethermin et al. (2010) prediction as a 
guide, our mm-selected sample never make up more than 
~ 50 per cent of the expected SFRd of IR luminous sources 
(Logio Lir > 12.5). We can see the reason for this clearly 
in Fig. [6] at Td = 30K we could expect our mm-sample to 
reach this level, while at Td = 40K we are only sensitive to 
IR sources above a luminosity of Logio Lir 12.7. Thus for 
the Bethermin et al. (2010) prediction to hold this shortfall 
must be made up of IR sources with Logio Lir <; 12.7 and 
T d > 30K. 


6 CONCLUSIONS 

We have presented an analysis of 63 radio-detected, mm- 
selected sources. Our parent sample of mm-sources was 
drawn from the recent MAMBO 1.2-mm imaging of Lock- 
man Hole North by Lindner et al. (2011) and the combined 
MAMBO-AzTEC 1.16-mm image of GOODS-N presented 
by Penner et al. (2011). Radio identifications are found for 
the Penner et al. (2011) sample, with counterparts for 24 of 
the 41 1.16-mm sources presented. These new identifications 
were combined with the 39 (out of 41) radio identifications 
presented by Lindner et al. (2011) for their MAMBO sample 
to give a sample of 63 radio-detected, mm-sources. We have 
cross-matched this mm-selected sample with SPIRE data 
from the HerMES project, making use of a map based tech- 
nique which reduces the problem of confusion in the SPIRE 
data by using radio and 24 /im positions as a prior. The main 
results of this work are as follows; 


Pir (logi 0 L Mpc" 3 ) 
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Figure 7. Contribution of mm-selected sources to the IR lumi- 
nosity density evolution. For comparison, we show the estimates 
of Chapman et al. (2005) and Wardlow et al. (2011) at z = 0— 
4. As ECDFS is underdense at submm wavelengths, Wardlow 
et al. (2011) present two estimates; one based on the observed 
number of sources in ECDFS, and one which has been scaled to 
the SHADES 850 /im number density (Coppin et al. 2006). Also 
shown is the predicted contribution by > 10 12 - 5 l© ULIRGs to 
the IR luminosity density from the empirical model of Bether- 
min et al. (2010). At z < 3 our estimate is significantly lower 
than both those found in both previous studies and models, while 
higher than previous studies at z ~ 4. These disagreements can 
be explained when considering the flux limits and biases of the 
comparison samples. 

• After comparing with the SPIRE data we found 47 mm- 
sources with a detection in one of more SPIRE band, giving 
a detection rate of 74 per cent. 

• For SPIRE detected mm-sources with spectroscopic 
or good quality optical/near-IR photometric redshifts, we 
found a tight correlation between dust temperature (Td) 
and IR luminosity (Lir). This correlation is found to be off- 
set to colder dust temperatures from the local relationship 
(e.g. Chapman et al. 2003). 

• Using the observed tight relationship between dust tem- 
perature (Td) and IR luminosity (Lir) to break the degener- 
acy between redshift and dust temperature we have demon- 
strated that photometric redshifts from submm/mm data 
alone offer an accuracy of |Az|/(l + z) = 0.16 (< |Az| >= 
0.51). If constraints from the radio-far IR correlation are in- 
cluded this is improved to |Az|/(l + z) = 0.15 (< |A«| >= 
0.45). 

• Via a combination of spectroscopic, optical/near-IR. 
photometric and submm/mm-radio photometric redshifts, 
the redshift distribution of mm-selected sources was pre- 
sented. We found that the redshift distribution of mm- 
selected sources peaks at higher redshifts than those selected 
at ~ 850 fi m, with a long tail to high -2 (z > 3.5) made 
up of primarily mm-only (i.e. SPIRE undetected) sources. 


This difference may be explained by a combination of the 
depth of the radio data used to identify the sources, as 
well as the evolving A-correction between submm and mm- 
wavelengths. Good agreement was found between this study 
and previous mm-selected samples. 

• We measured the contribution of mm-selected sources 
to the star formation rate density (SFRd) of the Universe. 
We found that mm-selected sources do not contribute sig- 
nificantly to the global SFRd, nor do they make up the bulk 
of IR luminous sources, at any redshift. 
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